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sured energy difference between the 2P and 2S states of muonic hydrogen disagrees significantly 
with that extracted from electronic hydrogen or elastic electron-proton scattering. This is the 
proton radius puzzle. The origins of the puzzle and the reasons for believing it to be very sig- 
nificant are explained. Various possible solutions of the puzzle are identified, and future work 
needed to resolve the puzzle is discussed. 
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1 Introduction 

The recent determination of the proton radius using the measurement of the Lamb 
shift in the muonic hydrogen atom [l] startled the physics world. The obtained 
value of 0.84184(67) fm differs by about 4% or 4.4 standard deviations from the 
electronic hydrogenate atom CODATA j2j value of 0.8758(77) fm and from a 
similar value with larger uncertainties determined by electron scattering [s] . The 
preceding sentence brings up a number of simple questions. The most prominent 
are: how can atomic physics be used to measure a fundamental property of a 
so-called elementary particle, why should muonic hydrogen be more sensitive to 
this quantity than electronic hydrogen, and why should a 4% difference between 
proton radii extracted using different techniques be important? In the present 
Introduction we sketch brief answers. The remainder of the article is devoted to 
detailed answers to these and related questions and implications. 

The sensitivity of atomic energy levels to the non-zero size of the proton is 
determined by the probability that the bound lepton be within the volume of 
the proton. This probability is roughly given by the ratio of proton to atomic 
volumes: {vp/aBY — (o^^rrp)^, where rp is the proton radius and rrir is the 
lepton reduced mass. The muon mass is about 200 times the electron mass so 
the muon is about 8 million more times likely to be inside the proton than the 
electron. 

To be a bit more precise we need to define the proton radius. The mean-square 
value Tp of the radius is given by 



2 _ (.^Ge 
rt. = — D 



where is the Sachs' electric form factor of the proton and is the negative 
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of the square of the four-momentum transfer to the proton. We note in passing 

that, because of recoil effects, this radius is not actually the integral of times 

a true density jijjsj . However it is that is observable in energy levels of muon 

and electron hydrogen atoms. 

For non-relativistic systems the lepton-proton Coulomb interaction is modified 

away from the point Coulomb interaction by the factor Ge{Q'^) which is the 

probability amplitude for the proton to absorb the exchanged photon (6|[7]. The 

resulting difference between the true potential and the potential for a point-like 

{Ge = 1) proton is given by 

SV(r) . VcM - V?(r) = -4™ / '^^'^J,'-" . (2) 

Using the expression G£;(q^) — 1 « — q^r^/G in Eq. (p| is a very accurate approx- 



imation because in atomic physics r-p g ~ rp/aB ~ 10^^. If this approximation is 
used, the perturbation 6V is a delta function at the origin multiplied by a factor 
that includes rp, and the resulting energy shift for atomic S-states is given by [S]: 

AE = {^s\SV\^s) = l7Ta\^s{0)fr^p. (3) 

The ^iid Pi/2 states are degenerate for solutions of both the Schroedinger 
and Dirac equations. But with the effects of non-zero proton extent included, the 
Si/2 state moves up and the P1/2 state is not affected. This gives a contribution 
to the energy difference between P1/2 states, part of the Lamb shift. 

The difference between the electronic and muonic determinations of the proton 
radius is very puzzling. One possibility is that experimental results are wrong. 
But multiple independent electron-proton experiments agree, and the muonic 
hydrogen experiment looks more convincing than any of the electron-proton ex- 
periments. Similarly, many measurements of transition frequencies in hydrogen 
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are in agreement. As detailed in Sect. 2.1, the muonic hydrogen charge radius 



would require a change of the Rydberg constant [2] by 4.9 standard deviations |T], 
in disagreement with the value obtained from many measurements in hydrogen. 
Another way out would be to assert that the various QED calculations that give 
contributions to the Lamb shift (other than the proton extent) are not suffi- 
ciently accurate. But all of the important contributions have been obtained by 
more than one group |6] and checked and rechecked. Another possibility is that 
the muon and electron have different interactions with the proton. But this vio- 
lates a well-known and well-tested principle called lepton universality. The failure 
to understand why the proton radius can be different in muonic and electronic 
hydrogen is indeed a puzzle. 

2 Electronic measurements 

Fig. [T] shows several determinations of the proton charge radius. Early elastic 
electron-proton scattering measurements from Orsay |9], Stanford [10] , Saska- 
toon [IT 12 and Mainz [l3], and the various re-analyses of these world data [l4| 
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. Spectroscopy of atomic hydrogen became sensitive to on the percent level 



in the mid-1990s. 

2.1 Hydrogen spectroscopy 

Spectroscopy of atomic hydrogen (H) and deuterium (D) has fueled the devel- 
opment of modern physics for more than a century. It was the discovery of 
the Lamb shift in H |17j which first showed effects of quantum electrodynamics 
(QED) which go beyond the Dirac equation. Today, the 1S-2S transition in H 



has been measured with an accuracy of 4 parts in 10^^ 18 . Other transitions 
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most notably the two-photon-transitions between the metastable 2S state and the 



states 8S,D 19 or 12D 20 , have been measured with accuracies around 1 part 
in 10^^. For a review of the relevant transition frequencies in H and D see [2]. 

QED describes the energy levels of H with extraordinary accuracy. The test of 
QED using measured transition frequencies in H is limited by two of the input 
parameters required in the QED calculations, namely the Rydberg constant i?oo 
and the rms proton charge radius rp. Hence, one can either supply any of these 
two numbers from a source other then H spectroscopy (such as rp from elastic 
e-p scattering of muonic hydrogen) and then test the correctness of QED, or use 
QED to extract the fundamental constants i?oo and rp. 

Somewhat simplified, the energies of S-states in H are given by 

E{nS)^-^ + '^ (4) 

where n is the principal quantum number, and Lis denotes the Lamb shift of 
the IS ground state which is given by QED and contains the effect of the proton 
charge radius rp. Numerically, Lis — (8172 -|- 1.56 rp) MHz when rp is expressed 
in fm, so the finite size effect on the IS level in H is about 1.3 MHz. 

The different n-dependence of the two terms in Eq. Q permits the determi- 
nation of both Roo and rp from at least two transition frequencies in H. Ideally, 
one uses the most accurately measured 1S-2S transition |18| and one of the 2S- 



8S,D/12D transitions 19,20 . The former contains the maximal IS Lamb shift 
Lis a-nd is hence maximally sensitive to rp. The latter contain only smaller 
Lamb shift contributions due to the 1/n^ scaling in Eq. Q , and hence determine 
Roo- Fig. [2] shows the difference values of rp obtained by combining the 1S-2S 
transition and each of the other precisely measured transitions in H. In addition, 
it contains three values of rp obtained from a direct measurement of the 2S-2P 
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transitions in H. These are not sensitive to the Rydberg constant. 

From Fig. [2] one can observe that all rp values from H favor a larger rp around 

0.88 fm. Still, half of the individual rp values agree with the muonic hydrogen 

value of 0.84 fm on the level of la. In fact, only the 2S-8D5/2 transition [19| 

disagrees with the muonic rp value on the level of 3a. The discrepancy between 

the combined value from H, as obtained in the elaborate CODATA adjustment 

of the fundamental constants [2], and the muonic hydrogen value, is about A.4a. 

2.2 Elastic electron-proton scattering 

Elastic electron scattering has been used to measure the electromagnetic structure 
of nucleons and nuclei for about six decades, and been reviewed for nearly as 
long [21[|22 . For the proton charge and magnetic radii, the key points are that the 



structure is encoded in two momentum-space, relativistic invariants - the electric 
and magnetic form factors GEiQ^) and GMiQ'^) ~ and the radii are determined 
from the slope of the form factors at four-momentum transfer = 0. (In this 
section the topic of the proton magnetic radius arises at times, but usually our 
remarks will be directed towards the charge radius.) In the one-photon exchange 
approximation, the measured experimental cross section is related to the form 
factors by 

da da 



dn exp dn point 



(5) 



1 + r 

Here r = Q'^/Am'^ and 9 is the lab scattering angle. Traditionally, cross section 
measurements were corrected for certain "radiative corrections" , higher order di- 
agrams beyond the one-photon exchange approximation, not included in Eq. [5j 
Form factors were extracted using Rosenbluth separations, multiple cross sec- 
tion measurements at different beam energies and scattering angles, resulting in 
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the same momentum transfer. Then fits of the form factors with simple, some- 
times theoretically inspired, parameterizations were used to determine the radii. 
However, due to a variety of issues with experiments and corrections, early radii 
extractions in the literature should no longer be taken seriously. One should in- 
stead be concerned with analyses of the past several years that pay attention to 
the issues mentioned below and are careful about extracting the radius. 

Many of the issues in extracting a proton radius from scattering data have been 
known for over a decade |14|[23| . These include 

• The treatment of systematic uncertainties in many experiments was often 
unclear, and apparently optimistic. Fits to world data, taking the uncer- 
tainties as reported, do not give reduced sufficiently close to 1. 

• The radius determination is mainly sensitive to low data, and world 
data fits can be satisfactory overall without going through the low data. 

• Fits extracting the radius using a single functional form typically ignored 
the uncertainty from model dependence. 

• In principle the Taylor series expansion Ge = 1 — Q'^r^/Q + Q^rp/120. . . 
gives a model independent radius determination. In practice the finite range 
of and uncertainties in the data lead to correlations between the coefficients 
Tp and possible truncation errors. Because the coefficients grow with 
order, there is no range where one term sufficiently dominates that its 
factor Tp can be fixed and used in fitting lower order terms to lower 
data. 

• Issues with the absolute normalization affect the extracted radius, but ffoat- 
ing the normalization decreases sensitivity to the radius, and increases sen- 
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sitivity to possible dependent errors - the rapid change of the cross 

section and data rates with is a problem. 

• Early radius extractions lacked Coulomb corrections, the acceleration of the 
electron in the proton's Coulomb field, which tend to increase the extracted 
radius about 0.01 fm 1231. 



The observation of differences in the ratio Ge{Q^)/Gm{Q'^) measured by Rosen 
bluth techniques and by polarization transfer techniques, first observed in [24] 
ultimately led to the recognition of the importance of two-photon exchange cor- 



rections 25 , ignored in earlier work. Calculations of two-photon exchange de- 
pend on the off-shell structure of the proton, and are not under precise theoretical 
control. Experimental measurements generally constrained two-photon effects to 
be no more than This topic continues to be under active investigation, due 

to its potential impact on knowledge of the proton structure, and of the radius. 
High precision experiments are under way to determine observables that depend 
on two-photon exchange, such as difference in e^p and e~p cross sections, single 
spin asymmetries, and variations in form factor ratios extracted at the same mo- 
mentum as a function of the beam energy or scattering angle. While two-photon 
exchange should be under better control in the near future, it is currently be- 
lieved to be a small correction for low data of most importance to the radius 
puzzle, and mainly of importance in determining the magnetic radius |26j . 

A renewed interest in the long-range, low-Q^ structure of the proton led to 
several new electron scattering experiments performed at the same time as the 



muonic hydrogen experiment. The Bates BLAST collaboration 27 used a polar- 
ized electron beam incident on a polarized gas target to determine that the form 
factor ratio I-lGe/Gm was close to unity for ~ 0.15 — )■ 0.6 GeV^. The Jef- 



10 Proton Radius Puzzle 

ferson Lab LEDEX collaboration [28| , |29j used polarization transfer to determine 

that the form factor ratio deviated significantly from unity by 0.3 - 0.4 GeV^, 
subsequently confirmed and improved upon in [30j. Note that while the polar- 
ization data give a form factor ratio, which does not determine the radius, the 
ratio does help constrain normalizations of cross section data during fits, leading 
to an improved value of the radius. The analysis of fSO] gives confirms the larger 
proton radius of electron measurements, Vp = 0.875 ± O.OOSexp i 0.006 ju- Data 
already taken to extend the form factor ratios to even lower Q^, 0.01 GeV^ - 0.06 
GeV^ are currently under analysis jsi] . 

Most importantly for the radius puzzle, the Mainz Al collaboration |^ mea- 
sured 1422 precise relative cross sections over a wide range of angles (20° - 135°) 
and beam energies (180 MeV - 855 MeV), corresponding to = 0.0038 0.98 
GeV^. The use of one of the three Al spectrometers at all times as a monitor 
spectrometer along with moving spectrometers a fraction of their angular ac- 
ceptance to new settings allowed high precision relative cross sections and form 
factors to be determined. The data were fit with a variety of functional forms 
using 31 normalization constants. Generally it was found that a good x^/d.o.f., 
~1.14, could be obtained with flexible forms (polynomials or splines) but not 
with inflexible forms, such as dipoles. Polynomials and splines led to slightly dif- 
ferent extracted radii, leading to rp = 0.879 ± 0.005stat =t O.OOAsyst ± 0.002modei 
± 0.004group) where the final uncertainty comes from the polynomial vs. spline 
difference. 

The fits of [3] and [30] use entirely independent data sets, and give form factor 
ratios in good agreement in the more limited range of [so] . The proton charge 
radii from the two fits also agree, but the magnetic radii are significantly different. 
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A large part of the difference appears to result from different treatments of the 



two-photon exchange corrections 32 , 33 , and it appears that the difference in 



magnetic radii has no effect on the determination of the electric radius. 

These new experimental reports reinforce the proton radius puzzle. Several re- 
examinations of issues in determining the radius by fitting scattering data have 
ensued. Here we highlight four recent contradictory results. 

Dispersion relations and the vector meson dominance model have long been 
used to fit the form factors, with the space-like form factors resulting from poles 
in the time-like region. Various other features such as continuum contributions 
and expected high-Q^ asymptotic behavior have been added. The most recent 



work of this sort 34 finds = 0.84 it 0.01 fm with x^/d.o.f ~ 2.2. Several notes 
are in order. Neutron and proton data are simultaneously fit. The Mainz Al 
data set normalization is allowed to float, and in most cases the change is less 
than 1%. Statistically, the reduced is inadequate; while in general this could 
result from underestimated uncertainties in the data, the lower of tbe flexible 
fits of [3,35 suggest this is not the case here. However, visually the fits of |34| 
go through all the lower data well, and start to diverge somewhat at larger 
values of . This suggests that for determining the radius the large value of 
is not an issue, but that there may be an issue with inappropriate strength in 
higher order terms - thinking in terms of the Taylor series expansion - leading 
to an incorrect result for the radius term. The smaller charge radius in these 



dispersion relation fits is in line with previous fits of this type 36 , 37 



The z expansion ~ see 38 for its application to the proton form factors - 
provides a method of incorporating physical constraints (analyticity) in choosing 
a functional form to fit the form factor data. The form factor is parameterized 
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as a power series in a complex variable z{Q^) that is constrained to lie within the 

unit circle. The power series expansion of the form factor is model independent, 

unlike the particular functional forms chosen in the dispersion relations analyses 



of 34 ,36 37 and others. The benefits of this transformation include a maximum 
value of z and coefficients that are bounded, guaranteeing a finite number of 
terms will be adequate. A fit with only the proton data yields rp = 0.870 fm it 
0.023 fm lb 0.012 fm; while including neutron data and the vrvr continuum leads 



to rp = 0.871 fm ± 0.009 fm ± 0.02 fm ± 0.02 fm - see [38] for details. While 
the analysis continues to suggest a larger proton radius, it also suggests that the 
proton radius uncertainties are underestimated by other recent analyses. 

Some of the issues discussed previously in relation to fits were revisited in recent 



work by Sick 39,40 . The extraction of the proton radius was found to be most 
sensitive to data in the range 0.01 GeV^ 0.06 GeV^ (0.5 fm~^ 1.3 fm"^). 
Fit functions with polynomials in a denominator can, with negative terms in the 
denominator, generate pathological behaviors - poles in the form factors at high 
outside the range of data. Thinking in terms of the Fourier transform of a 
spatial distribution, this generates non-physical oscillatory behavior that extends 
to large r and affects the radius extracted. Despite the low of the Mainz data, 
fits of the data would be more reliable if the data extended to lower and / 
or did not have a floating normalization. Parameterizations that correspond to a 
more sensible large r fall off should be safer. A sum-of-Gaussians r-space fit gave 
stable radii, without the issues discussed, of rp = 0.886 fm it 0.008 fm. 

Recent unpublished work by C.E. Carlson and K. Griffioen points out that a 
linear fit to the lowest Mainz data, < 0.02 GeV^, where Ge appears to be 
entirely linear, yields rp ~ 0.84 fm, rather than ~ 0.88 fm. But great care must 
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be taken in doing any such fit 14 . Since for the proton > 0, the quadratic 

term in the Taylor series expansion nearly guarantees that a linear fit to low 

data will underestimate the radius. An estimate of the potential size of the effect 



can be made with the Kelly form factor parameterization 41 , which has ~ 
0.86 fm. A fit to pseudodata up to 0.02 GeV^, with the density and statistics of 
the Mainz data set, yields ~ 0.84 fm it 0.01 fm. Adding in a quadratic term 
leads to uncertainties on rp too large to distinguish between 0.84 fm and 0.88 
fm. Similar concerns apply to an extension of this approach up to 0.2 GeV^. An 
extensive study of this issue has been performed by M.O. Distler. 

To summarize, the apparently simple problem of determining the slope of the 
form factor at = has numerous potential pitfalls. The weight of the evidence 
at this point continues to favor a larger radius, about 0.88 fm, but suggests that 
claiming an uncertainty at the 0.01 fm level is optimistic. 

3 Measurement in muonic hydrogen 

A measurement of the Lamb shift in muonic hydrogen (//p) was initially con- 



sidered half a century ago as a test of electron vacuum polarization effects 42 
complementary to the Lamb shift in regular, electronic hydrogen [l^ which is 
dominated by the electron self-energy 143] . The first observation of x-rays from 



muonic hydrogen succeeded shortly afterwards 44 . 

Laser spectroscopy of the 2S-2P transition in /ip requires, however, |Up atoms in 
the metastable 2S state. Several groups failed to observe such long-lived /ip(2S) 
atoms when muons are stopped in molecular hydrogen gas [45| - |48| , see also [49| . 
The first observation of long-lived fip atoms in the 2S state |50| - [52) was the 
starting point [53 - 55 of the recent Lamb shift measurement [T 56 . 
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When negative muons are stopped in molecular H2 gas at low pressure (1 mbar 

at room temperature), about 1% of the muons form /ip(2S) atoms |50[|57] with 

a lifetime of about 1 /xs [52] . At higher gas pressures the 2S state is quickly 

de-excited in collisions with H2 molecules. 

A novel beam line for low-energy negative muons |55| was built at Paul- 

Scherrer-Institute (PSI) in Switzerland which delivers about 600 /i' per second 

with a kinetic energy between 3keV and 6keV. About half of the muons stop 

in 1 mbar H2 gas within a target vessel which has a length of 20 cm along the 

muon beam axis. The transverse dimensions of the muon beam are 0.5 x 1.5 cm^. 

Hence the muon stop volume is small enough to be illuminated with laser light 

of sufficiently high fluence (see below). The muons arrive at random times, so 

each muon has to be detected before it enters the target vessel. Two stacks of 

ultra-thin carbon foils, separated by ~ 35 cm, are used to detect the muon with 

high efficiency (80% and 70%, respectively) using secondary electrons ejected by 

muons crossing the carbon foils. A signal in both carbon foil stacks, with the 

correct time-of-flight for a low-energy muon traveling the distance between the 

carbon stacks, serves as a start signal for the data acquisition and the pulsed 

laser system. 



The pulsed laser system 158 59 delivers 5 ns long pulses of light with a wave- 



length tunable from 5.5 ^um to 6.1 /im, with a pulse energy of about 0.25 mJ. A 



cw-pumped Yb:YAG disk laser system 59 , 60 produces 100 mJ of pulsed pump 



light only a few hundred nanoseconds after a muon-induced trigger signal at a 
random time. This light is used to amplify red light, tunable around 708 nm, 
using Ti:sapphire crystals. The red pulses are converted to the desired infrared 
(IR) wavelength at 5.5-6 nm via three sequential Stokes shifts in a high-pressure 
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Raman cell |61j. The frequency of the laser light is stabilized and controlled 

at 708 nm, and calibrated in the IR using well-known absorption lines of water 

vapor (H2O). These H2O lines are known with an accuracy of a few MHz, but 

pulse-to-pulse energy fluctuations of the laser system cause a broadening of the 

laser line width resulting in the final 300 MHz uncertainty of the laser frequency 

calibration [I]. 

A non-resonant multi-pass mirror cavity [62] inside the target gas vessel ensures 
efficient illumination of the muon stop volume inside the H2 target gas. The 
estimated laser fluence is about 6mJ/cm^, which results in approximately 30% 
transition probability when the laser is tuned to the center of the resonance. 

The successful 2S-2P laser excitation is signaled by the emission of a 1.9 keV 
Kq x-ray which is emitted in the radiative 2P-1S de-excitation that follows the 
2S-2P transition immediately. The x-rays are detected in 20 large-area avalanche 
photo diodes (LAAPDs, 14 x 14 mm^ active area each) 63,64 which are mounted 
close to the muon stop volume. 

The experiment [l] records time spectra of Kq, x-rays for various laser frequen- 
cies (i.e. wavelengths) [65]. For each time spectrum, the number of Kq, x-rays 
in the time window at which the laser illuminates the muon stop volume is nor- 
malized by the "prompt" x-rays which are emitted by the 99% of the muons 
which do not populate the metastable 2S state, but proceed directly to the IS 
ground state. 

The resonance curve in Fig. [3] is obtained by plotting the normalized number 
of laser-induced x-rays as a function of laser frequency. The fit of this reso- 
nance curve reveals the resonance frequency and hence the Lamb shift in muonic 
hydrogen [l]. 
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4 Potential solutions 

The various possible explanations of the proton radius puzzle are discussed in 
this section. 

4.1 Muonic hydrogen experiment 

The discrepancy between the observed resonance position in muonic hydrogen [I] 



and the predictions using the Vp values from elastic electron scattering 14,16] or 
hydrogen spectroscopy |15j (which together form the CODATA values 2006 and 
2010) is enormous, as can be seen in Fig. [3] 

The resonance center is 75 GHz away from the central value of the CODATA 
prediction. This corresponds to about four times the width of the resonance. The 
statistical uncertainty of the resonance center is about 700 MHz, i.e. 1/100 of the 
discrepancy. The systematic uncertainty of the measurement is 300 MHz. It is 
due to the laser frequency calibration mentioned in Sect. [3j Other systematic 
effects like ac- and dc-Stark shifts, Zeeman shifts, Doppler shifts etc. have been 
investigated but were found to be even smaller [I]. 

It is obvious from these numbers that the discrepancy is severe, and it is difficult 
to imagine an effect that could shift the resonance position by 4 line widths. 
Jentschura [66] found that the presence of an electron could result in a shift of 
the resonance position if the distance between the electron and the /ip(2S) atom 
was about 1 (electronic) Bohr radius. He suggested that the spectroscopy might 
have happened not on a ;up(2S) atom, but on the molecular ion pfie~ . 

Karr and Hilico [67] studied the 3-body systems pfj,e~ and pp^~ and found that 
these are not stable. In addition, the molecular ion ppfj,~ has been calculated 
to be very short-lived |68|-l70|. The observation of the quenching of long-lived 
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/ip(2S) atoms |52| supports this picture. Hence, the proton radius puzzle cannot 

be explained by 3-body physics and molecular ion formation. 



4.2 Theory of muonic hydrogen 

Muonic hydrogen, an atom consisting of the proton and the muon is very similar 
to regular "electronic" hydrogen, so the theoretical description of the two atoms 
has many elements in common. The main difference is in the mass of the lepton. 
The muon is about 207 times heavier than the electron, and the muon-proton 
mass ratio r], according to NIST [2], is given by r/ = 0.1126095272(28). This is not 
very small as is the case for eH, and therefore the static nucleus approximation 
is not necessarily the best starting point for the theoretical description of ^H. 
Nevertheless, the traditional nonperturbative approach introduced by Borie and 
Rinker in a seminal work [tT], continued by Borie |72[|73| , and very recently by 
Indelicato [74] and Caroll et al. [75] is based on the Dirac equation and corrected 
for the finite nuclear mass by including some additional contributions. Here, we 
present an approach, symmetric in the proton and the muon, that can also be 
applied to electronic hydrogen. There are no significant numerical differences 
between both approaches [7^, but the one presented below, to our opinion, is 
simpler and more elegant, since all the corrections can be accounted for in a 
systematic way. 

The binding energy ii^ of a two-body system, according to quantum electrody- 
namics (QED) is a function of the fine structure constant a and of the muon- 
proton mass ratio r] = m^j^/mp: 



where we assume that both of the constituent particles are point-like. We neglect 




(6) 
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also, but only for a moment, the presence of electron-positron pairs which intro- 
duces another dimensionless parameter. Even after these simplifications, there 
is no single equation, which gives the exact binding energy, as is the case for 
the Dirac or Schrodinger equation. To some extent, the Bethe-Salpeter equation 
can be regarded as the most general bound state equation, but the kernel of this 
integro-differential equation can be given only perturbatively in the fine structure 
constant, and thus the exact solution is not known. In the approach presented 
here, traditionally called nonrelativistic quantum electrodynamics (NRQED), one 
expands the binding energy £(a,rj) in powers of the fine structure constant a, 
with = 137.035999074(44) and derives an exact closed formula for the ex- 
pansion coefficients. In most cases these coefficients are expectation values of 
some effective Hamiltonian, and thus can easily be calculated. 

The description of the QED bound two-body system begins with the leading 
order part, which is the nonrelativistic energy of two particles interacting via the 
Coulomb force in the reference frame with vanishing total momentum: 

z I nip r 

According to the Schrodinger equation the binding energy is a function of the 
principal quantum number only 

mrO^ , . 

^o(n) = (8) 

where rur stands for the reduced mass, and fi = c = \. If so, the difference 
A£J = Els in binding energies of 251/2 and 2Pi/2 states 

Els = E2P,^, - = 0. (9) 

Departures of Els from zero, traditionally called the Lamb shift, therefore arise 
from additional corrections to the energy. 
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The first set of corrections we discuss involves relativistic effects, although these 

are not the largest ones. In the static nucleus approximation the binding energy 
of 2Si/2 aiid 2Pi/2 states according to the Dirac are also exactly the same. More- 
over, the leading recoil corrections also do cancel, therefore relativistic effects in 
muonic hydrogen ~ ry^ and thus are quite small. How are these relativistic cor- 
rections calculated? According to the NRQED approach, the leading relativistic 
corrections can be expressed in terms of the so called Breit-Pauli Hamiltonian [77| 

Hbp = Hq + SHbp 

^»BP = -^-fj- P- f- + ^) P- (10) 
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where s^, Sp are spin operators of the muon and the proton, and g^, gp = 
5.585 694 712(46) are ^-factors. The interaction 6Hbp includes also the effect 
due to the finite proton charge radius rp together with the Darwin terms in order 
to make it clear that these two effects are separated. Hbp was derived years ago 
before NRQED was formulated, nevertheless we use it as a part of the NRQED 
approach. The related correction to energy, neglecting terms with the proton 
spin Sp, the proton charge radius rp and assuming g^ = 2^ (corrections due to 
gfj_ — 2 are included later as the muon self-energy) 

5,,iEls = {2Py2\SHBp\2Py2) - {2Sy2\SHBp\2Sy2) 



A S 

a m 
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0.05747 meV (11) 



is quite small, as expected. This means that the perturbative treatment of rela- 
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tivistic effects is very appropriate and higher order terms wih be even smaher, if 



not negUgible. In deriving (11), we neglected the hyperfine interaction. Although 
we are not interested here in the hyperfine splitting, this interaction leads to the 
mixing of Pif2 ^3/2 states and shifts these levels. This additional mixing 



correction was first considered in 78 and amounts to 



5ELs{P[/f) = -0.U46meY, (12) 

and ELsiP[/2^) = +0-1446 meV. By definition, we do not include this as part of 
the Lamb shift, but consider it as a separate shift of Pj^^ levels, see Fig. |4[ 

So far we showed that the relativistic effects are relatively small. Indeed, the 
leading effect on the Lamb shift comes from vacuum polarization effects. Namely, 
the Coulomb interaction between the muon (electron) and the proton is modified 
by creation of the electron-positron pairs in the electric field. In the language 
of quantum electrodynamics the photon propagator includes closed fermion loop 
corrections, what results in the modification of the Coulomb interaction by the 
so called Uehling potential 



Kp(r) = -^« f^^e— n(,^). (13) 



where 

The characteristic range of the electron vacuum polarization Vvp{r) is about the 
electron Compton wavelength, which is close to the Bohr radius in ^uH. Namely 
the ratio /3 = rrif./ {rrir oi) = 0.737386 is not far from one. This means that the 
parameter (3 must be kept intact in making expansions in a, because the muon 
nonrelativistic wave function has a significant overlap with the electronic vacuum 
polarization potential V^p. As a result, the vacuum polarization correction is large 
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for muonic hydrogen, in contrast to that in the electronic hydrogen. Indeed, the 

leading vacuum polarization correction 

S.pELs = (2Pi/2|Kp|2Pi/2) - (25i/2|Kp|25i/2) = 205.0073 meV (15) 

is the dominating part of the muonic hydrogen Lamb shift f78|, and all other 
corrections are at least two orders of magnitude smaller. Note, that the expec- 
tation value is taken with nonrelativistic wave function and the muon-proton 
mass ratio r/ is included exactly. At first glance this procedure is less accurate 
than the expectation value taken with Dirac wave function, but we point out 
that the use of the reduced mass in the Dirac equation is more than question- 
able, because the leading relativistic recoil correction is not properly accounted 
for by using the reduced mass. In the present discussion, the vacuum polariza- 
tion effects are treated perturbatively by taking expectation values instead of 
solving numerically the Schrodinger equation with this potential. Higher-order 
effects obtained by solving the Dirac equation can be included in the perturbative 
approach by including the higher-order terms as well. The slight advantage of 
perturbative approach is that the vacuum polarization potential acquires higher 
order QED corrections (two- and three-loops), therefore it is also known pertur- 
batively, so higher order perturbation and loop corrections are treated on equal 
footing. Namely the second order correction due to the Uehling potential, 

SEls = 0.1509 meV, (16) 

is combined with the two-loop vacuum polarization (one-particle reducible and 
irreducible two-loop diagrams) 

5^L5 = 1.5081 meV. (17) 

The last correction is clearly the dominant one. Although much smaller, the 
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complete correction due to electronic vacuum polarizations is included also in 

present //H theory. It was calculated, first by Kinoshita in |79|, and later slightly 

corrected by Karshenboim et al. in Ref. 180]. The final result is 



= 0.0053 meV. (18) 

The higher order (4-loops or more) vacuum polarization corrections are negligible 
at present. There are however additional vacuum polarization corrections, which 
are quite interesting. The electron-positron pairs are distorted by the presence of 
the real particles. Namely, these pairs can additionally interact with the muon or 
with the proton. This leads to three box type of diagrams, and the overall effect 



as obtained by Karshenboim et al 81 



5Els = -0.00089(2) meV, (19) 

is almost negligible, as these three diagrams tend to cancel out. We nevertheless 
include them in the overall theoretical predictions to demonstrate that vacuum- 
polarization effects are calculated as completely as possible, and no significant 
effects are being neglected. 

There are also vacuum polarization correction from muon pairs, and this effect 
is included later together with the muon self-energy. Finally the hadronic vacuum 



polarization, considered in detail in the context of the muon g — 2, amounts to 82 



5hvp^L5 = 0.0112(4) meV. (20) 

Up to this point, all the vacuum polarization corrections has been calculated 
using the nonrelativistic wave function. The relativistic O(a^) effects combined 
with evp (electron vacuum polarization) are included separately as another cor- 
rection. Again, the muon-proton mass ratio r/ is treated exactly, so the use of 
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Dirac equation is not appropriate. Instead, we introduce vacuum polarization 

effects in the Breit-Pauli Hamiltonian. This is achieved as follows. One notes 



from Eq. (13) that the evp corrected Coulomb interaction can be represented 
as a Coulomb interaction obtained with the exchange of massive photons, inte- 
grated with the mass dependent weight function. So, one derives the modified 
Breit-Pauli Hamiltonian d^pHsp obtained using the massive photon, neglects 
the hyperfine interaction and finds that the corresponding correction to binding 
energy is 

5wp,vc\ElS = i^vpHBp) + 2 (Kp TT^ JTTT -f^Bp)- (21) 

[1^ - N) 

The resulting contribution to the Lamb shift was first calculated in [78] , however 



with some errors. The correct result was first obtained by Jentschura in 83 and 
confirmed by Karshenboim in |84| and amounts to 

(^vp.rei^^LS = 0.01876 meV. (22) 

If one used the Dirac equation in the infinite nuclear mass limit, the obtained 



result 0.021 meV would differ significantly from ( 22 ) . Nevertheless, this correction 
is quite small, as any relativistic correction for //H, so higher order relativistic 
effects combined with evp can safely be neglected. 

What next? In electronic hydrogen the dominating part of the Lamb shift is 
the electron self-energy. In the case of /xH it is a small but nevertheless important 
correction. The formula for the one- loop Lamb shift in eH is well known |6], the 
resulting correction in fiH is obtained by replacement of the electron mass by the 
muon mass with the result 

(5E = -0.6677 meV. (23) 
Further related corrections arise from including both muon self-energy and elec- 
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tron vacuum polarization effects. The calculation is not very simple as it involves 

modifying Bethe logarithms by evp. The complete calculation was performed by 

Jentschura [85] and the resulting correction is 

= -0.0025 meV. (24) 

We pass now, to remaining corrections which have an overlap with the proton 
elastic structure and polarizability effects, and not always are treated consistently 
in the literature. The pure recoil corrections of order a^, derived originally by 
Salpeter, assume that both particles have spin 1/2 and are point-like. They can 
not be obtained from the Dirac equation by including the reduced mass, as the 
Dirac energy is even in a. Their derivation requires a full QED treatment and 
the obtained result [g] 

17( A "^r jZaf r 2 / 1 \ 8 1 7 



5/0 



where 



^n/ V 2 ■'■ n/ 2nJ ' 1{1 + 1){21 + 1) 

is valid for an arbitrary mass of orbiting particles. Later, in the next section, 
we consider the elastic and inelastic two-photon exchange contribution, and this 
Salpeter correction has to be consistently subtracted out from the elastic ampli- 
tude. Coming back to pure recoil effects in //H, the resulting correction is 

SEls = -0-04:50 meV. (27) 

The last correction we consider is due to the proton self-energy and the related 
definition of the proton charge radius 1781. We mention in passing that for an 
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arbitrary spin / nucleus the charge radius is defined in NRQED by the effective 

coupling 

where E is the electric field external to the proton. For a spinless nucleus = 0, 
for half-spin 61/2 = 1/8, what corresponds to the so called Darwin- Foldy term. We 
point out for future considerations that for higher spins there is no unique value 
of 61. Here, the muon is point-like, so its charge radius by definition vanishes, 
and the charge radius of the proton, the main issue of this review, is defined by 
the above equation. This definition would be equivalent to that obtained from 
the slope of Ge{Q'^) form factor at = 0, if the electromagnetic form factors 
are meaningful at an arbitrarily high precision level. This is not necessarily the 
case because of QED radiative corrections. The proton self energy leads to the 
modification of elastic form factors in such a way that they depend on a fictitious 
photon mass, or in other words, are no longer well defined. The only correct 
definition has to employ complete two-photon structure functions, but this will 
not be analyzed here. We take the simplest possible point of view and use the 
formula for the low energy part of the proton self-energy 

The corresponding correction in /xH is 

5Els = -0.0099 meV . (30) 

The high energy part of the Lamb shift is by definition included in the pro- 
ton formfactors, more precisely in the charge radius and the magnetic moment 
anomaly. The problem which remains is how this definition corresponds to the 
proton charge radius as obtained from the electron scattering. We do not aim to 
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analyze it, because possible inconsistencies are much smaller than the observed 

discrepancy between fiH and the electron scattering determination of r^. 

All corrections presented in this section sum to 

A£;ls = 206.0330 meV. (31) 

There are many further small QED corrections considered in the literature which 
are higher order in a or in powers of vacuum polarization. The main purpose of 
these further calculations was to find a possible missing large correction. How- 
ever, no significant QED correction has been found yet, which would explain the 
discrepancy, but in principle we can not exclude its existence. 

All known corrections have been summarized recently by Antognini et al. |76|, 
and we refer readers to that work for more details. Their value (Eq.(|32[) below) 



differs only slightly from that in Eq. (31). 



The calculation of the fine and hyperfine splittings [86], see Fig. |4] except for 
IS 

-^HFS^ is much simpler, so we present here only the final numerical results 



76 



A£;ls = 206.0336(15) - 5.2275(10) + A^tpe 

A£^FS = 8.3521 meV 

A£;Jpg/' = 22.8089(51) meV, (exp. value) 

A£;hfs = 7.9644 meV 

= 3.3926 meV 

A = 0.1446 meV (32) 

where A£^tpe = 0.0351(20) meV (see Table 1) is a proton structure dependent 
two-photon exchange contribution, which is considered in the next section. The 
theoretical value for 25^/2 hyperfine splitting is uncertain due to not well known 
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proton-structure effects in magnetic interactions, so the above value is the exper- 
imental one, taken from 56 . 



4.3 Proton Polarizability Effect in Muonic Hydrogen 

The previous subsection is concerned with the QED calculations needed to deter- 
mine the proton radius from the muonic energy level splittings between 2P and 
25" states. There is one particular term, occurring at the intersection between 
QED and strong interaction effects which is deserving of special attention. This 
is the proton polarizability contribution that enters in the two-photon exchange 
term, see Fig. [5] . The computed effect of this term is proportional to the lepton 
mass to the fourth power, and so is capable of being relevant for muonic atoms, 
but irrelevant for electronic atoms. 

The proton can be excited by the absorption of a photon and then de-excited 
by the emission of a photon. Thus the two-photon exchange term depends on the 
forward virtual photon-proton (Compton) scattering amplitude T^^'^{v^q^) where 

is the square of the four momentum, of the virtual photon and v is its time 
component. The contribution of the graph of Fig. [s] is given by the expression 87 



e 



4 



d^q 1 



E = --m'J ^^[r^'^-t'^^(M)]V(m), (33) 
where t^'^ is the Compton amplitude for a pointlike fermion of mass M (proton) 

3 3 

or m (muon), 0^(0) = "g™'' for the 2S state with m,. as the lepton reduced mass. 



The role of the term t^'^{M) in Eq. (33) is to subtract the term recoil correction 



discussed in the previous section. Using Eq. (33) still requires the subtraction of 
the infrared divergent term related to the QED correction to the proton charge 
radius discussed in the previous section. 

The quantity T^'^{v, q^) is decomposed into a combination of two independent 
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terms, Ti^2{i^,q'^), allowed by symmetries. The imaginary parts of the functions 

Ti 2{y,q^) are related to structure functions Fi^2 measured in electron- or muon- 
proton scattering, so that Ti^2 can be expressed in terms of F12 through dispersion 
relations. However, Q^) falls off too slowly for large values of u for the dis- 
persion relation to converge. Hence, one makes a subtraction at = 0, requiring 
that a new function of Q^, Ti(0,Q'^), be introduced. The contributions to the 
Lamb shift are broken into three terms [88 1: E = AE'^'- -|- Ai?*"^' -|- Ai?*"**, where 
el refers to lepton-proton intermediate states, inel refers to intermediate involving 
excitations of the proton, and subt refers to the subtraction term. This separa- 
tion into three terms is somewhat arbitrary, as emphasized by Ref. [89| . But so 
far no authors have found that differences between different possible separations 
are relevant for resolving the proton radius puzzle. 

The subtraction term has typically been handled by using a power series ex- 
pansion around = 0, and then using effective field theory to determine the 



coefficients of the series. In practice, people 87,88,90,91 have used 

Ti(0,Q2) = ^q2^ (q2)^ . (34) 

Q 

where -Fioop(Q^) is a function that falls with sufficiently rapidly so that the 
necessary integral converges. The calculations of three groups are summarized 
in Table [T] The theoretical uncertainty as determined from the three values of 
AE shown in the Table would seem to be very small. However, there is a hint of 
trouble in the large relative differences between the values for Ai?''"^*. Birse & 



McGovern 91 reevaluated the subtraction term (using the separation of terms 
employed by Ref. [88]) and computing one order higher in the power series for 
to find AE"''''* = 4.2±1.0//eV. Some authors |89||92(]94| have argued that the 
uncertainty in this term is substantial, even though the typical values are only a 
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few /ueV, very small compared to the 310 fi eV needed to account for the proton 

radius puzzle. 

The trouble is that the integral that determines A£^*"^* j|94j would diverge 
logarithmically without the introduction of the form factor -Fioop(Q^)- In princi- 
ple, one needs to evaluate an infinite set of terms in the power series describing 
Ti(0, Q^) and hope that the result would lead to a convergent integral. Miller | |94| 
used the form 

The use of this form allows one to state the expression for the energy shift in 
closed form, which is well approximated by 



A£--«-ma</,^(0)^. (36) 

Using, = 0.5 GeV and a = 0.0078/(4m^) =0.177 GeV^^ led to value of the 
Lamb shift of 0.31 meV. A relatively large value of 1/a must appear to obtain 
the necessary contribution, but no fundamental principles have been violated. 

Thus the above discussion shows that procedures used to estimate the size of 
the subtraction term suffer from significant uncertainties. This arises because the 
chiral EFT is being applied to the virtual-photon nucleon scattering amplitude, 
and as a result the energy shift depends on an integral that would be logarith- 
mically divergent if the function -Fioop were replaced by unity. Another way to 
proceed would be us to use an effective field theory (EFT) for the lepton-proton 
interaction. In EFT logarithmic divergences identified through dimensional regu- 
larization are renormalized away by including a lepton-proton contact interaction 
in the Lagrangian.This was done |94| by using standard dimensional regulariza- 
tion (DR) techniques by evaluating the scattering amplitude of Fig. [5j The term 



of interest is obtained by including only Ti{0,Q ) of Eq. (34) with Fioop = 1- 
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The result is : 



^(loop) = -i a^m— [- + log ^ + - 7^; + log A-K]ufUiU fUi, (37) 



-t a m - + log ^ + - 

2 a e b 



where lower case spinors represent leptons of mass m, and upper case proton of 
mass M, is Euler's constant, 0.577216- • • , /x is the parameter introduced in 
using dimensional regular ization and one works in a space of dimension d = A — e. 



The result Eq. (37) corresponds to an infinite contribution to the Lamb shift in 
the limit that e goes to zero. In EFT one removes the divergent piece by adding a 
lepton-proton contact interaction to the Lagrangian that removes the divergence, 
replacing it by an unknown finite part. The finite part is obtained by fitting to 
a relevant piece of data. Here the only relevant data is the 0.31 meV needed to 
account for the proton radius puzzle. When this is done the result is 

47r 



MJ^^ = i 3.95 a^m^UfUiUfUi. (38) 

X 



where the results are expressed in terms of the chiral symmetry breaking scaleA^ = 
4:11 Jtt, (/tt is the pion decay constant). The coefficient 3.95 is of natural size. Thus 
standard EFT techniques result in an effective lepton-proton interaction of nat- 
ural size that is proportional to the lepton mass. 



The present results, Eq. (36) and Eq. (38) represent an assumption that there 



is a lepton-proton interaction of standard-model origin, caused by the high- 
momentum behavior of the virtual scattering amplitude, that is sufficiently large 
to account for the proton radius puzzle. This assumption needs further testing. 
Fortunately, our hypothesis can be tested in an upcoming low-energy fi^p, e^p 
scattering experiment that is discussed in Sect. 5.2. 
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4.4 Theory of hydrogen energy levels 

One of the possible explanations for the proton charge radius puzzle is that it is 
caused by a mistake or a missing correction in the hydrogen Lamb shift theory. 
We will argue that this is the least probable solution. Why? It is because 
corrections to the hydrogen Lamb shift have been calculated by many groups 
using different methods, with only few exceptions which are described below in 
detail. For a recent review of Hydrogen theory in the context of determination 
of the proton charge radius see Ref. [2]. 

All the corrections to the Lamb shift are classified in powers of a, Z a and 
the mass ratio mg/mp, where the power of a counts the number of loops. The 
additional proton structure corrections, elastic and inelastic two-photon contri- 
butions are of the order of tens of Hz for IS state [2], while the discrepancy in 
the proton charge radius corresponds to about 100 kHz. So, any proton structure 
effect beyond the the finite size correction is negligible. 

The leading contributions to the hydrogen Lamb shift, the one- loop electron 
self-energy and the vacuum polarization have been calculated by an expansion 



in Z a 95 or directly numerically 96 , using the known form of the relativistic 
electron propagator in the Coulomb field. Both methods are in agreement, but 
the numerical one is presently more accurate, about a few Hz for the ground 
state, so there is no room for any mistake of order of 100 kHz. 

The two-loop contribution has also been calculated numerically f97| and ana- 



lytically 95 , 98 , but the results differ slightly. Numerical calculations have been 
performed only for the nuclear charge Z > 10 and the result for Z = 1 was ob- 
tained by extrapolation. The analytic calculation was numerically accurate, but 
some higher order terms were omitted, because they are too difficult to evaluate. 
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The final result is taken as an arithmetic mean with an uncertainty that matches 

both values. This uncertainty corresponds to about 2 kHz for IS state, much too 
small to explain the discrepancy. 

The three loop contribution in the leading order of Za has been calculated 
only by one group |99|. This correction can be represented in terms of three-loop 
electromagnetic form factors of the free electron. Although technically very ad- 
vanced, the conceptually simple methods of the evaluation of Feynman diagrams 
can be applied, so we do not expect a mistake here. The binding corrections, 



next order in Z a are known only partially 100 , but can be roughly estimated 
to be about ~ 1 kHz. 

The pure recoil effects, two-body diagrams without self-energy and vacuum 
polarization, although small, are not very simple to evaluate because the full 
QED formalism must be used. The derivation of the closed formula for the 



pure recoil correction was achieved independently by two groups 101 102 . The 



numerical evaluation was performed by analytic expansion in Z a 102 , 103 and 



directly numerically 104 . Both are in agreement. 

The radiative recoil corrections have been calculated only by Za expansion. 



but independently by two groups 105 , 106 . The overall result is rather small, 
some logarithmic higher order terms are also known j[2j, and remaining higher 
corrections are irrelevant at the order of discrepancy. 

In conclusion there is no room for any new correction, or any mistake in the 
present calculations that would shift IS Lamb shift by about 100 kHz. 
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4.5 Physics Beyond the Standard Model 

The possibility that the proton radius puzzle might be caused by a difference 
between the fxp and ep interactions is very exciting because such an effect is not 
part of the Standard Model (SM). Awareness of the possibility of such a failure of 
what is known as electron-muon (or more generally lepton) universality was re- 



awakened by the muon g — 2 experiment 107 -109 . The experimental value of the 



muon anomalous magnetic moment exceeds the SM expectation by more than 
three standard deviations: Aa^ = aJfP — aj;*^ = (29 it 9) x 10^^°. This difference 
could be caused by a new interaction. In addition, new gauge forces mediated 
by particles in the MeV-GeV scale could be dark matter candidates [110 -112 



Obtaining a theory or model that could account for the proton radius puzzle, the 
muon anomalous magnetic moment and predict other new verified phenomena 
would be a tremendous breakthrough. However, it is not easy to invent new 
interactions that differentiate between electrons and muons because the principle 
of lepton universality is very well tested. 

The most relevant constraints of lepton universality have been nicely summa- 
rized in Refs. [113 114 . The basic ideas are to use the non-observation of new 



particles that couple only to muons and to search for differences between muonic 
and electronic decays of various particles. Constraints are obtained from the de- 
cay of the T resonances [115 -117 ; neutron interactions with nuclei 118 - 120| ; the 



anomalous magnetic moment of the muon 109 ; x-ray transitions in Mg and 



^^Si atoms; J/^ decay 121 ; neutral pion decay 122 ; and eta decay 123 



Barger et al. [113 postulate the existence of new interactions mediated by the 
exchange of a boson between muons and nucleons. To explain the proton radius 
puzzle, these exchanges must give rise to an attractive interaction that accounts 
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for 0.31 meV. One can include the exchange of scalar, vector and tensor bosons. 

Axial vector exchange gives a spin-spin interaction that affects the hyperfine 

splitting, but not the Lamb shift, and is not considered further in Ref. |113j . 

Given the non-relativistic motion of the muon and the nucleon, the exchange 

of boson between the muon and nucleon leads to a Yukawa interaction : 

AVir) = -ax^^> (39) 

where is the mass of the postulated boson and is the product of lepton- 
and nucleon-boson coupling constants divided by Air. Assuming sufficiently weak 
couplings, first order perturbation theory can be used to obtain the contribution 



to the muon Lamb shift 71 



as 



6{AE) = a^m^^^^Yf^' ^^^^ 



arrir ' 



where the muon- nucleon reduced mass is m^. 

Barger et al. conclude that new spin 0,1 and 2 particles that mediate flavor- 
conserving non universal spin-independent interactions are excluded. This finding 
depends on assuming that the new particles couple only to the muon, and that 
the coupling to mesons is the same as to the nucleon. 



Brax &; Burrage 124 note that light scalar mesons coupled to matter appear 
in many fundamental contexts, such as in the inflationary scalar field and in 
attempts to unify gravity with the standard model. However, these interactions 
are strongly constrained by experimental searches for fifth forces, and violations 
of the equivalence principle. Non-linear dynamics [125 -127 can be introduced to 



avoid these constraints. Brax & Burrage use the assumption that matter couples 
universally to the scalar field to obtain a new stringent bound on the coupling of 
scalar fields to matter from the precise measurements of the difference between 
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the energies of the Is and 2s states in hydrogen, while the couphng to photons is 

constrained by the Lamb shift. The net result of these considerations is that the 

contribution of a scalar field to the discrepancy between proton radii measured 

in electronic and muonic atoms is negligible. 

Tucker-Smith & Yavin |128 go through a chain of reasoning similar to that 

of Barger et al. fll3| , and find that a new force carrier with mass of order MeV 

that couples to protons and muons can account for the proton radius puzzle as 

well as the difference between QED and measurements for g — 2 oi the muon. 

They acknowledge that many constraints exist, but state that none of the existing 

constraints exclude such a force in a model-independent manner. The difference 

in conclusion with Barger et al. arises because Tucker-Smith & Yavin do not 

assume universal couplings to nucleons and mesons. Tucker-Smith & Yavin state 

that their postulated force can be tested in measurements of muonic-deuterium 

and helium. 

Motivated by the interest in in the possibility of new gauge forces the medi- 
ate interactions between dark matter and standard model parameters, Batell et 



al. 129| use a model in which a new U{1) gauge boson (heavy photon) is kineti- 
cally mixed with hypercharge. New scalar and vector force carriers interact with 
right-handed neutrinos. Current experiments at JLab, Mainz, KLOE and BaBar 
are searching for such particles but haven't found any. This limits the value of the 
mass of the heavy photon to be about 30 MeV. If the strength of this particular 
new force is large enough to account for the proton radius puzzle, parity violat- 
ing asymmetries in muon-nucleon or muon-nucleus scattering would be enhanced 



by several orders of magnitude and would be of order 10 . Ref. 114 use the 



nonobservation of missing mass events in the lepton kaon decay K — )■ /ix to find 
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strong constraints on parity- violating gauge interactions of a right handed muon 

provided x decays invisibly or does not decay inside the detector. 

Carlson &: Rislow [130 consider two models, one containing scalar and pseudo 



scalar particles and the other vector and axial vector particles. The basic idea 
is that the scalar or vector particles account for the Lamb shift while the pseu- 
doscalar and axial vector particles do not contribute. Large effects in the muon 
g — 2 are prevented by cancellations between the effects of scalar (vector) and 
pseudoscalar (axial) mesons. Constraints on the masses of the new particles are 
obtained from the nonobservation of such particles in kaon decays. Exchanges 
of pseudoscalar and axial vector particles contribute to hyperfine splittings, and 
the resulting constraints are not considered in Ref. [130| . 

It is also worthwhile to discuss the possible influence of new forces on the PSI 



muon-proton scattering experiment (Sect. 5.2). If the new particles are light, 
they act like a very small modification of the Coulomb force and can not be 
detected. If the new particles are sufficiently heavy, the force between the muon 



and the proton is well-approximated by a contact interaction (as in the model 94 
discussed in Sect. 3e) and its effects are detectable as an interference between a 
short-ranged and long-ranged Coulomb interaction. 

To summarize, there are many constraints on new forces. Some of these con- 
straints are derived from model-dependent considerations, and the new force 
possibility remains viable for a subspace of all of the relevant parameters. The 
ultimate disposition of this possibility will depend on either directly detecting 
the postulated particles or directly verifying an independent prediction of one of 
these models. 
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5 New projects 



37 



There is an interesting set of new experiments that could shed hght on the proton 
radius puzzle. 

5.1 New electron scattering experiments 

Improving the precision of the radius determination from electron scattering ex- 
periments requires that cross sections be measured with higher precision and/or 
to smaller Q^. All higher precision measurements that have been done to date 
have used magnetic spectrometers with small or moderate solid angle acceptance, 
and it is hard to see how a significant improvement can be made in the precision 
of these measurements. As ~ AE E' sin^ {9 / 2) , going to smaller requires 
smaller scattering angles and/or beam energies. But the magnetic spectrometers 
have limited angular ranges, and the experiments already run have used the low- 
est beam energies available at modern accelerators. Thus, an improved electron 
scattering experiment requires a different technique. 

The only experiment currently approved to improve the electron scattering 



radius determination is Jefferson Lab E12-11-106 131 , which aims to extend the 
range of electron scattering experiments from the 0.0038 GeV^ of the Mainz 
experiment down to about 10~^ GeV^. The experiment has several novel features 
that combine to allow a significant radius determination from a measurement 
only at low Q^. Low current, 1 and 2 GeV beams in Hall B impinge upon an 
open-ended hydrogen gas cell, as has been used in storage rings previously, to 
eliminate the need to subtract target end cap backgrounds. The spectrometer, 
the PRIMEX HYCAL (PbW04) calorimeter, will cover the entire angle range 
of ~10 mr - 70 mr in one setting, to cover the range 10"'^ GeV^ — )• 10~^ 
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GeV^. MoUer scattering of the beam from atomic electrons will provide the data 

used to determine the detector solid angles, offsets, and efficiencies at the high 

level needed. The lower energy of the MoUer scattered electrons allows the two 

reactions to be distinguished by the calorimeter. 

The difference in the cross sections arising from the difference in the measured 
radii is small, requiring high precision cross section measurements. Approximat- 
ing Ge = ^ — radii of 0.842 fm vs. 0.875 fm lead to differences in the 
cross sections of 0.005%, 0.05%, 0.5% and 1.0% for of lO"'' GeV^, lO'^ GeV^, 
0.01 GeV^ and 0.02 GeV^, respectively, and differences half that large for the 
form factor. Existing electron experiments each are about 5a different from the 
muonic hydrogen result. If the radius is to be determined to give a similar qual- 
ity result, then relative cross sections need to be determined at the 0.2% level. 
The large cross sections and event rates for low scattering are large, making 
statistical uncertainties at the 0.1% quite feasible. 

Systematic uncertainties must also be controlled at a level that would be very 
difficult without a well understood calibration reaction. Reaching low requires 
small scattering angles and/or beam energies. With a 1 GeV beam energy, reach- 
ing 10~^ (10~^) GeV^ requires a scattering angle of 10 (32) mr. But the Mott 
cross section diverges at small angles due to the 1/ sin^(^/2) factor, leading to 
a shift in the cross section of 0.01%//xr (0.03%//xr) at 10"^ (lO""^) GeV^. A 10 
/xr knowledge of the scattering angle is then needed to limit shifts the relative 
cross sections by 0.2%. In contrast, typical high precision electron scattering ex- 
periments with magnetic spectrometers determine absolute scattering angles at 
about the 1 mr level. 

Jefferson Lab E12-11-106 is approved and tentatively scheduled to run in 2014- 
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2015. If various technical challenges are met, it has an excellent chance of re- 
solving all experimental questions related to data normalizations and fits used to 
extract the proton radius from electron scattering experiments, either confirming 
the puzzle, or resolving it as being due to incorrect radius extractions from the 
electron (scattering) data. 

5.2 Elastic muon-proton scattering 

The missing measurement in determining the proton radius with muons and 
electrons is muon scattering, a measurement that the MUon proton Scattering 
Experiment (MUSE) collaboration proposes at PSI [132 . Muon scattering is 



not an easy measurement, as muons are produced as secondary beams from the 
decays of vr's produced from a primary proton beam. As a result, muon beams 
have low flux and large emittance, compared to electron beams. Unless efficient 
particle separators are available, the muon beam will also include background 
pions and electrons. 

The proposed MUSE experiment is based at the PSI vrMI beam line. The 
beam line provides of order MHz muon fluxes for momenta from about 100 —)• 
250 MeV/c. To handle the multiple particle species in the beam, MUSE uses 
Scintillating Fiber arrays to determine RF time and only operates at momenta 
where the muon RF time is about 4 or more ns different from electron and pion 
RF times. A custom FPGA system will allow the RF time to be determined in 
hardware so that the information is available to the trigger. To avoid issues arising 
from accidental coincidences, MUSE intends to limit the channel acceptance to 
keep the total beam flux at ~5 MHz. The large emittance of the beam requires 
the use of GEM chambers to track individual beam particles into the target. 
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so that the scattering angle is sufficiently well known. Test measurements have 

largely verified that the beam properties are sufficient for the experiment. 

The low rates also require a large acceptance spectrometer. MUSE is planning 

a nonmagnetic spectrometer, consisting of wire chambers and fast scintillators for 

triggering and particle identification. The main reason for this is that experiments 

with large acceptance detectors and magnetic fields do not measure precise cross 

sections. 

MUSE intends to measure with both and beams, so that two-photon 
effects can be directly determined from the data. MUSE also intends to measure 
both muon and electron scattering at the same time, so that the experiment can 
compare them. 

The collaboration has studied singles rates, background triggers, and system- 
atic uncertainties to estimate how well the cross sections, form factors, and radii 
from the various measurements can be compared. A number of difficulties arise 
that are not present in electron scattering experiments. The tt induced back- 
grounds are efficiently removed at the trigger level through the measurements of 
beam RF time, which is only possible due to the low beam flux. The beam muons 
decay, which is mainly a problem as the decay electrons cannot be distinguished 
from the muon scattering at the trigger level, and are not as easily distinguished 
with a non-magnetic spectrometer. Although the muon decay fraction is small, 
it is much larger than the scattering fraction. The muon radiative tail is very 
small, but electrons have a large radiative tail which must be handled. The beam 
emittance is large, requiring measurement by detectors in the beam, which fur- 
ther increases the emittance. Many of the systematic uncertainties are similar 
to electron scattering experiments, but the low energy muon beam has reduced 
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sensitivity to angle offsets, radiative corrections, and hadronic backgrounds com- 
pared to a GeV electron beam. Preliminary estimates are that the radius can 
be extracted from muon scattering about as well as in the Mainz experiment, in 
part because the experiment can reach ~ 0.002 GeV^, about half the lower 
limit of the Mainz measurement. 

The MUSE experiment requires significant new funding (~2 M$) to be carried 
out. It might be able to run as soon as 2016. 

5.3 Spectroscopy of electronic atoms and ions 



As detailed in Sect. |2.1l the value of Vp can be determined from the 1S-2S transi- 



tion in H [1^ using precisely calculated QED theory only if a sufficiently precise 
value of the Rydberg constant Roo is supplied. At present, Roo originates from 
several measurements of transition frequencies in H and D. 

In fact, the correlation coefficient between R^o and rp is 0.984 in the 2010 
CODATA adjustment [2j. Hence, any new determination of the Rydberg constant 
at a level of a few parts in 10^^ would help to shed new light on the proton radius 
discrepancy between electronic and muonic hydrogen. 

Several such projects are underway. Flowers et al. at the British NPL are 
measuring the 2S-6S/D 2-photon transitions in atomic H [133 . Beyer et al. at 



MPQ, Garching, aim at an improved value of Roo from 1-photon transitions in 
H, namely the 2S-4P and higher P-states. This project utilizes for the first time 
a cryogenic beam of H atoms from a nozzle at 6K, optically excited to the 2S 
state using the 1S-2S laser system |18| . This can eliminate possible systematic ef- 
fects like dc-Stark shifts from patch charges caused by the electron bombardment 
traditionally used to create H(2S) atoms and align the H atom beam onto the 
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laser beams. The cold atom beam can help to reduce systematic effects caused 

by the large velocity of traditional "hot" beams, but the 1-photon transition is 

particularly sensitive to 1st order Doppler shifts. The 1S-3S transition in H has 

recently been measured by Nez et al. at LKB, Paris [134 . It is the second best 

measured transition frequency in H, but as can be seen in Fig. [2] it does not 

provide the most precise rp value when compared with 1S-2S following Eq. Q. 

This is because both the 1S-2S and the 1S-3S transition depend on the large IS 

Lamb shift Lis, and the "lever arm" on R^o is hence reduced. Nevertheless, this 

measurement is very valuable because it serves as a cross-check for the 1S-2S 

transition in H ^18j which has for a long time been measured by only the MPQ 

group. The Paris group is improving their 1S-3S setup, in particular the laser 



system, and aims at a significant improvement of their accuracy 134 . The 1S-3S 
transition in H is also being measured at MPQ, Garching, by Peters and cowork- 



ers 



135 . To overcome difficulties in the generation of the 205 nm light required 



to drive the 1S-3S transition in H they use two counter-propagating picosecond 
frequency combs. 

Laser spectroscopy of He atoms and H-like He"*" ions are pursued by Udem et 



al. at MPQ 136 and by Ubachs, Eikema and coworkers in Amsterdam |137|. 
Ultimately, this will allow a Rydberg constant determination using He atoms or 
ions. A precise value of the He nuclear charge radius is available from electron 
scattering experiments 



138 . Soon, the Lamb shift in muonic helium ions 139 



will improve this accuracy further (see Sect. 5.4). 

An alternative route towards a new determination of i?oo is pursued by Tan 
and coworkers at NIST 140 . Here, H-like highly charged ions (e.g. Ne^''") are 
created in circular large-£-states at high principal quantum number n and stored 
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in a Penning trap. In large-£-states the effect of the finite nuclear size on the 

energy levels is negligible, and measuring these highly charged ions at high n 

states ensures that the transition frequencies can be determined using lasers in 

the optical region. 

Hessels et al. at York University, Toronto, aim at a new measurement of 
the "classical" 2S-2P Lamb shift in H. This radio-frequency transition was last 



measured in 1994 141 . The new direct measurement of the Lamb shift will yield 
a value of rp that does not require the knowledge of i?oo , so it is an independent 
cross-check of the 2S-ni transitions in H and D that determine Roo today. 

5.4 Spectroscopy of exotic atoms 

Spectroscopy of exotic atoms will also contribute to the solution of the pro- 
ton radius puzzle. Precision laser spectroscopy of purely leptonic systems like 

can test 
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positronium (Ps = e~^e~) 142,143 or muonium (Mu = /^^e~) 
bound-state QED free from finite nuclear size effects. Improved spectroscopy of 
the 1S-2S transitions and the ground state hyperfine splittings are on the way 
or being prepared. Ps and Mu test QED in the electronic and muonic sectors, 
respectively. 

Laser spectroscopy of muonic deuterium atoms will soon shed new light on the 



radius puzzle 145 . A new charge radius of the deuteron can be compared with 



the one obtained from the isotope shift of the 1S-2S transition of electronic H and 



D 146 . The muonic deuterium measurement may also give new insights into the 



question of nuclear polarizability (see Sect. 4.3). 

The PSI project R-10.01 is going to measure the Lamb shift in muonic helium 
ions /i^He and //^He 139 . The measurement will give a tenfold improved value 
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of the charge radii of the stable hehum isotopes, compared to the value from 



electron scattering 138 . Together with the new measurements of transitions in 



electronic He atoms or ions 137 the muonic helium Lamb shift will provide a 
cross-check of the proton radius puzzle at Z=2. 

A future possibility may be the measurement of transitions in muonic Li, Be 



and B ions 147| . Accurate charge radius differences of these nuclei are known 



from electronic isotope shift measurements [148 -150 



6 Conclusions 

The proton radius puzzle is real. A physical parameter should not depend on 
the method of extraction. Yet here we have two highly precise methods based on 
electronic and muonic hydrogen atoms that disagree significantly. The possible 
origins of the disagreement are the subject of much of this text. To summarize: 
the possibilities are 

• The electronic hydrogen experiments are almost, but not quite, as accurate 
as stated (Sect. 2.2). 

• The QED calculations are almost, but not quite, as accurate as stated (Sect. 
4.2). 

• The two photon exchange term that depends on proton polarizability has 
not been correctly evaluated (Sect. 4.3). 

• The electron and muon really do have different interactions with the proton 
(Sect. 4.4), so that there is physics beyond the Standard Model. 

None of these possibilities seem very likely, but all must be pursued. 
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There are some tasks for theorists. While the results of many of the QED 

calculations have been reconfirmed, it is still possible that something is missing 

in the theory. Current theories of beyond the standard model physics are very 

primitive and need technical improvements before being realized as a complete 

gauge theory without anomalies. Better treatments of the two-photon exchange 

term could be very useful. 

Much additional work by experimentalists is necessary, and any future progress 
is contingent upon continuing a strong experimental effort in diverse directions. 
The electronic hydrogen experiments need to be redone using modern techniques 
(Sect. 5.3). A new high energy-low momentum transfer electron scattering mea- 
surement could reach the accuracy of the spectroscopic determinations and pro- 
vide an independent confirmation or repudiation of the radius obtained from 
electronic hydrogen spectroscopy, Sect. 5.1. The proposed muon-proton scatter- 
ing experiment (Sect. 5.2) would rule on theories of the two-photon exchange 
interaction and also on theories beyond the standard model. Discovery of the 
new particles required by theories beyond the standard model would make such 
theories believable. Exotic atoms, in which muons are bound to nuclei, would 
also rule on all possible theories. Sect. 5.4. Future work on positronium would 
put QED to detailed tests. 

The proton is the only stable baryon and all of Life depends on it. It is 
important to understand the proton. Its radius should be a simple quantity to 
determine and understand, but this is not the case. Recent experimental results 
have not been understood. But future work will either reveal the true value of 
this radius, a better understanding of its structure, or a very unexpected feature 
of its interactions. 
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Figure 1: Proton radius determinations over time. Electronic measurements seem 
to settle around rp=0.88fm, whereas the muonic hydrogen value |1| is at 0.84 fm. 



Values are (from left to right): Orsay M, Stanford 10, Saskatoon 11,12 



Mainz \13\ (all in blue) are early electron scattering measurements. The red sym- 
bols between 1996 and 2001 originate from laser spectroscopy of atomic hydrogen 
and advances in hydrogen QED theory (see [2j and references therein) . The green 
and red points in the year 2003 denote the reanalysis of the world electron scat- 
tering data [I4] and the world data from hydrogen and deuterium spectroscopy 
which have determined the value of rp in the CODATA adjustments [2}[ 
the 2002 edition. 
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Figure 2: Proton charge radii rp obtained from hydrogen spectroscopy. According 
to Eq. Q, Tp can best be extracted from a combination of the 1S-2S transition 



frequency [18| and one of the 2S-8S,D or 12D transitions [19 20 . The value from 
muonic hydrogen |l] is shown with its error bar. 



Table 1: Numerical results for the 0(a^)m^ proton structure corrections to the 
Lamb shift in muonic hydrogen. Energies are in yueV. 

(/ueV) 



Ref 



Ref. 



78 



87 



Ref. 



90 





5.3 ± 1.9 


1.8 


2.3 




-12.7 ±0.5 


-13.9 


-16.1 




-29.5 ± 1.3 


-23.0 


-23.0 


AE 


-36.9 ±2.4 


-35.1 


-36.8 
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Figure 3: Resonance in muonic hydrogen, together with the positions predicted 
using the proton radii from elastic electron-proton scattering using pre-2009 world 
data 14,16 and the CODATA-2006 value from H spectroscopy 15. 
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Figure 4: Level scheme of muonic hydrogen for n = 2 shell, the artificial 2P]^/2 
and 2P3/2 levels corresponds to centroid for A = 0. Numerical values for level 



splittings is presented in Eq. ( 32 ) . 




Figure 5: The box diagram for the ©(a^m^) corrections. The graph in which the 
photons cross is also included. 



